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I. INTRODUCTION 



A key property of the statistics of quantum systems in two space dimensions is pro- 
vided by the possibihty to display intermediate fractional statistics interpolating between 
bosons and fermions: the properties of anyons, the two-dimensional identical particles obey- 
ing fractional braiding statistics and carrying fractional charge, have been the subject of 
an intense and continuing interest jl-^. Both Abelian and non-Abelian anyons (associated 
to respectively one- dimensional and irreducible higher-dimensional representations of the 
braid group) have been extensively studied both for their intrinsic interest and their con- 
nection with quantum Hall systems j5-7|. In particular there is an increasing interest in the 
investigation of the properties of non-Abelian anyons for their application to topologically 
fault-tolerant quantum information processing jsj. 

The ideal gas of anyons is also interesting for the phenomenon of statistical transmuta- 
tion {2-^, i.e. the fact that one can treat non-interacting anyons as interacting bosons or 
fermions. The idea behind statistical transmutation is that one can alternatively consider 
the system as made of interacting particles with canonical statistics or of non-interacting 
particles, but obeying non-canonical statistics. This makes in general difficult the study of 
the equilibrium thermodynamical properties of the ideal anyon gases j4| and therefore any 
result that gives even qualitative informations on such ideal gases is valuable. 

Dating back to the seminal work by Uhlenbeck and Gropper 8|, a standard way to 
characterize the effects of the quantum statistics on the properties of ideal gases is provided 
by the determination of the so-called statistical potential. As detailed in textbooks 
one can show that the partition function (PF) of a gas of particles approaches - for sufficiently 
high temperatures - the PF of the classical gas (with the correct Boltzmann counting). When 
this computation is done for an ideal quantum gas (under the condition that the thermal 
wavelength is much smaller than the interparticle distance) one finds that the quantum PF 
becomes the PF of the classical ideal gas. Evaluating the first quantum correction, one 
can appreciate that the quantum PF can be formally written as the PF a classical gas in 
which a fictitious two-body interaction term (the statistical potential) is added {^[[lol. The 
statistical potential gives a simple characterization of the effects of the quantum statistics 
of the ideal gases: the statistical potential is "attractive" for bosons and "repulsive" for 
fermions, and it is respectively monotonically increasing (decreasing) for bosons (fermions). 
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■ the statistical interparticle potential 



io|. 



Another important result is that a suitable integra]_of 
gives the second coefficient of the virial expansion 

The statistical potential of a gas of ideal Abelian anyons has been studied in |lll-ll3| and 
it depends on the statistical parameter a (we remind that a = and a = 1 corresponds re- 
spectively to free two-dimensional spinless bosons and fermions, while a = 1/2 corresponds 
to semions jj]). It is found that for 1/2 < a < 1 the statistical potential Va{r) is mono- 
tonically decreasing, while for 0<a<l/2it has a a minimum at a finite value of r and 
it is increasing for larger values of r (while for a = is monotonically increasing). We can 
refer to these behaviors respectively as "quasi-fermionic" (1/2 < a < 1) and "quasi-bosonic" 
(0 < a < 1/2). The purpose of the present paper is to compute the statistical potential of 
the ideal gas of non-Abelian anyons: we find that the behaviour of the statistical potential 
depends on the Chern-Simons coupling and the isospin quantum number. As a function of 
these two parameters, quasi-bosonic and quasi-fermionic regions emerge, and they are part 
of a phase diagram which will be presented below. 

The plan of the paper is the following: in Section 11 we recall the steps which lead 
to the computation of the statistics potential Va of an ideal gas of Abelian anyons with 
statistical parameter a. Since it is possible to show that the statistical potential for the 
non-Abelian gas may be written in terms of sums of Abelian statistical potentials (having 
different statistical parameters depending on the projection of the isospin quantum number), 
we provide in Section II a detailed study of Va- we present a compact and useful integral 
representation for it, showing that it can be written in terms of bivariate Lommel functions 

i n 

14( 1 . Furthermore, using the Sumudu transform of the statistical potential [15[, we also 
give a closed expression of Va{r) in terms of the inverse Laplace transform of an algebraic 
function of r and a. Using this result we are able to give a simple expression for the 
statistical potential of the ideal gas of semions, and we show that for a general value a it is 
possible to retrieve the well-known result for the second virial coefficient of an ideal anyon 
gas found in 1^ (with an hard-core boundary condition for the two-body wavefunction at 
zero distance). For the sake of comparison with the non-Abelian case that follows, the limit 
behaviors both at small and large distance of the statistics potentials are presented. In 
Section III we introduce the non-Abelian Chern-Simons (NACS) model studied in the rest 
of the paper: we compute the statistics potential (within the hard-core boundary condition 
frame) as a function of the Chern-Simons coupling k and the isospin quantum number I 
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and we build a phase diagram summarizing the behavior of the statistical potential in terms 



of K and /. We then show that the second virial coefficient, previously studied in |17H20|. 
is correctly retrieved. The asymptotic expressions for the small and large distance of the 
statistics potentials are also given. Finally, our conclusions are drawn in Section IV, while 
supplementary material is presented in the Appendices. 

II. STATISTICAL POTENTIAL FOR ABELIAN ANYONS 

In this Section we introduce the model for an ideal gas of Abelian anyons, and we then 

as a function of the statistical parameter a, obtaining 



derive its statistical potential Va{r\ 
the expression for Va{r) given in llMlSj], and also providing an explicit formula for the 
semions (half-integer values of a). The results for Va{r) and the asymptotic expressions for 
small and large distance will be used in the next Section, where the statistical inter-particle 



potential of an ideal gas of non-Abelian anyons is derived and studied. 



fl. 



Abelian anyons admit a concrete representation by the fiux-charge composite model 
and the statistics of these objects can be understood in terms of Aharonov-Bohm type 



interference 
anyons reads 



21 



22| . The Hamiltonian for the quantum dynamics of an ideal system of 



n=l 

where p„ is the momentum of the ra-th particle (n = 1, ■ ■ ■ , A^). Similarly we will denote 
the position of the ra-th particle by r„ = In Eq.(Il]) a is the statistical parameter: 

notice that the physical quantities, e.g. the virial coefficients, are periodic with period 2 j4|: 
the bosonic points are defined by a = 2j and the fermionic ones by a = 2j -|- 1, j integer. 
For this reason we will consider in the following a G [0, 2]. 

In Eq.([T]) a„ is the vector potential carrying the flux attached to the bosons: indeed, the 
Hamiltonian ([1]) is written in the so-called bosonic representation, i.e. it is an Hamiltonian for 
identical bosons, and therefore acting on the subspace of wavefunctions which are symmetric 
with respect to the exchange of particles. The explicit expression for = (a^, a^) is 

I ' n. ' 



where e*-' is the totally antisymmetric tensor = 1,2). 
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Let's recall that this model of Abelian anyons also admits a field-theoretic description: 
in fact (non-relativistic) anyons can be described by bosonic Schrodinger fields ip and 



coupled to a Chern-Simons gauge field living in (2+l)-D [23|, l2j| (then = 0, 1, 2). The 
Lagrangian density of such a system reads 



where c gives the measure of the interaction among particles mediated by the U{\) gauge 
potential a^, with the covariant derivatives given by Dt = dt + iqao, D = V — iqa, and the 
anyon statistical parameter a to be identified as a = g^/(27rc). 

In the study of a quantum-mechanical ideal gas, the effect of the symmetry properties 
of the wave function can be interpreted, from a classical point of view, as a consequence 



a 



referred to as 



of a fictitious classical potential introduced by Uhlenbeck and Gropper 
effective statistical potential, which represents the first quantum correction for the classical 
PF jsl, [lo| . For our purposes we have to consider the two-body case, which is relevant for the 
subsequent computation of the statistical interparticle potential p, [loi . The statistical po- 
tential completely determines the second virial coefficient, which gives the thermodynamical 
properties of the system in the dilute (high-temperature) regime. For the two-anyon system, 
after separating in ([1]) with N = 2 the center-of-mass dynamics (i.e. that of a particle having 
mass 2M), one is left with the dynamics of the relative wavefunction: 

1 i2 
- ^ [V,. - z a a(f)J ip{r} = E ij{r} , (3) 

where r = ri—r2 is the relative coordinate, and a*(r) = e*-'^. Therefore the relative dynamics 
is equivalent to that of a single particle in presence of a point vortex a(r) placed at the origin. 
An analysis of the statistical interparticle potential for an ideal gas of Abelian anyons is 



presented in [ll|, The eigenf unctions of reads 

where the capital (italic) letters respectively refer to center-of-mass (relative) coordinates. 
The bosonic description used in ([T]) imposes the condition I = even; furthermore Jm{x) 
denote the Bessel functions of the first kind [l^ (their definition is recalled in Appendix Rll . 
Notice that the wavefunctions (jlj) are the eigenfunctions of ([3]) provided that the hard-core 
boundary conditions \l/a(0) = are imposed. 
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The two-body PF is Z = Tre~^^^ = 2AX^'^Z', where Z' is the single-particle PF in the 
relative coordinates, (3 = l/ksT, At = (/3/i^/27rM)^/^ is the thermal wavelength and A is 
the area of the system. The relative PF Z' is given by 



(5) 



l=—oo 



(with p = Hk). It is possible to conveniently rewrite Eq.(|5]) by using the following integral 



property 



25(1 of the Bessel functions: 



1 



e-"V, (2/3v^) J, (27v^) dx = - L 



-(/32+72)/a 



(6) 



iQ a \ a J 

where ly is the modified Bessel function of the first kind (see also Appendix |A]). The 
relative PF then takes the form 



Z' 



-. oo 



dxe ""/i 



\2n-a\[X) , 



where 



X 



2f3f9 ^ XI 



(7) 



(8) 



The one-body PF for classical systems, used in (8|] to define the concept of effective statistical 
potential, is 

^ ' (9) 



therefore, comparing ([7]) and (|9]) produces as a result [11 1 

oo 
n=—oo 

(see Appendix [B] for details). Eq. flTII]) is plotted in FigJT] 

For the considerations which follow, it is convenient to introduce the function 

oo 

Ma{x) = ^ I\2n-a\{x) , 

n.=— oo 

SO that the statistical potential can be written as 



(10) 



(11) 



(12) 



The inter-particle statistical potential admits a closed expression in terms of the bivari- 
ate Lommel functions (alias, Lommel functions of two variables). Indeed, as evident from 
Appendix lAl it is 



M-cXx) = i " Ua{ix, ix) — i" U2-a{ix, ix) 



(13) 
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Exp[-;Sv(r)] 
2.0 




0.1 0.2 0.3 0.4 0.5 0.6 

FIG. 1. Plot of e~^^"(^^ vs. t/Xt for different values of the statistical parameter: from top to the 
bottom it is Q = 0, 0.1, 0.2, 0.4, 0.5, 0.6, 0.7, 1. 

where Ua denote the Lommel functions of two variables [l^. Notice that the symmetry 
property Va{r) = V2-a{r) (Wa G TZ ) holds for the statistical potential. 

In Appendix [B] we prove the following integral representation for Va{r), which will result 
useful in Subsection III Al for discussing the large- distance limit behaviour: 



-/3i)«(r) 



-2x 



COS an 



sm air 



TX 



sinh [(1 - 

sinh t ^ ' 



Other integral representations for e"^""^*"^ are presented in Appendix [Bl In the remaining 
parts of this Section we discuss the asymptotic behaviour of for small and large values of 
the dimensionless parameter x, and we give a closed expression for the statistical potential 
Va_[f) in terms of the inverse Laplace transform of an algebraic function of r and a: this 
manipulation allows us to straightforwardly regain the known expressions for Vair) in the 
bosonic/fermionic cases, to obtain its expression in the case of semions and to finally recover 
the value of the second virial coefficient for a generic a presented in 16|. 



A. Limit behaviours 

In order to quantitatively understand the tendency of anyons to bunch together or vice 
versa to repel each other in different limit regimes of density, let's recall and further discuss 
ihe asymptotic behaviors of their effective statistical potential, for small and large distances 



12 



13|. 



8 



For small r (that is x ^ 1) we can approximate the summation term in (fTO!) as 



|2n-a| \^) 



n=0 



(^/2) 



+ 



r(2n + a + 1) ' r(2n + 3- 
[r(a + \)Y\xl2Y + [r(3 - a)]-i(x/2)2-- 



a) 



(x/2) 



2n+2-« 



Since 

/3^;,(r) = - 

it follows 



In 



2 5^ /, 



|2n.-a|ia^j 



In [2 e-"([r(a + l)]-i(x/2r + [r(3 



(15) 

a)]-i(x/2)^-")] 
(16) 



In [2 - 27rr V A|] - In 2 + f (r /At) 



(17) 



« = 0, 2 

ln[2 (7rrV2A^)"/r(a + l)], 0<a<l 
In [27rrVA^] , a = 1 

^ - In [2 (7rrV2 A2,)2-"/r(3 - a)] , 1< a < 2 

We may summarize the small distance behaviour as follows: Vaij) is repulsive and log- 
arithmically divergent to +oo for any a G (0,2), whereas for a = 0, 2 it is attractive, and 
quadratically increasing in r starting from the finite value v^ir = 0) = — ln2. The small-r 
asymptotic function for fa(r) is discontinuous in a = 1, being twice than the limit asymp- 
totic functions for a ^ \^ (in fact two equal Bessel terms contribute to the asymptotic 
behaviour for a = 1, whereas only one of them dominates when a 

To study the behavior of the statistical potential for large distance r (i.e. x ^ 1) 
we employ the integral representation (fT^ of the modified Bessel function. The method 
of steepest descent allows to evaluate to an arbitrary order the last term of this integral 
representation. At the first significant order, we get 

1[oL — 1) sin avr 



smavr 



TT 



sinh [(1 - a)t\ ,„,ht ^ 



-2z 



sinh t V27ri 
Therefore the large distance behavior of the statistical potential is given by 

V2(l 



(18) 



cos ttTT + 



a sm avr 



Trr/A] 



-27rrVA|, 



(19) 



for any a G [0,2]. Let us notice that this result differs from the corresponding one in 12], 
and that the asymptotic behaviours for a = 0, |, 1 (see formulas (1261) . ( 127|) . (!28|) in the 
sequel) are correctly retrieved. The statistical potential for large distance is vanishing for 
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ySv,(r) 
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FIG. 2. f3va[r) vs. r/Xx for different values of a: from top to bottom it is a = 
1, 0.7, 0.5, 0.3, 0.2, 0.1, 0. The bosonic curve (a = 0) is monotonically increasing wliile the fermionic 
curve (a = 1) is monotonically decreasing and divergent for r — t- 0. All the curves for < a < 1/2 
diverge for r — t- and have a minimum point. 

r — oo, and the interval a G [0, 1] (as well as the interval a G [1,2], due to the symmetry 
property Va = f2-a) is divided in two regions: for large distance, Va{r) is attractive for 
< a < 1/2, and repulsive for 1/2 < a < 1. 

The large-distance and short-distance behaviors, considered together, imply that Va{r) 
must admit a minimum point at finite distance, Vcria) for any < a < 1/2 (see Fig 12]). We 
denote the corresponding dimensionless quantity by Xcr(tt) = vrr^j,(a)/A^. The minimum 
point Xcr(a) tends to -|-oo for a — | , as shown in FigOl 

Figj2] clearly shows that Abelian anyons have, from the point of view of the statistical 
potential, a "quasi-bosonic" behavior for < a < 1/2 (i.e., Va{r) is non-monotonic with 
a minimum) and a "quasi- fermionic" behavior for 1/2 < a < 1 {va{r) is monotonically 
decreasing without a minimum). 

B. Laplace transform of the statistical potential and the 2"'^-virial coefficient 

In this Section we write an explicit formula for e~^'"°'^^^ as the inverse Laplace transform 
of a function of x and a. This result will allow us to write a simple formula for the statistical 
potential for the semionic gas and to compute the second virial coefficient, which of course 
coincides with the result reported in the seminal reference fl6| . 



Let's start by writing down the Sumudu transform [15|| of the function A^q, defined in 
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-2.5 -2.0 



' logf- 
-1.5 -1.0 V2 



a 

2 f 



FIG. 3. The dimensionless value of the minimum point Xcr(ct) of the statistical potential of 
Abelian anyons as a function of the statistical parameter a near a = 1/2. 

Eg. (ITT!) . The Sumudu transform of M.^ is defined as 

/■OO /"OO 7i /"OO 

\M^{x)\s = / e^'^Ma{0x)d9 = / e-*/"A^„(t)- = s / e-^*M,(t)dt , (20) 
Jo Jo ^ Jo 



where s = -, whence 



Ma - 

X 



xC[Ma]ix) 



(21) 



J s 



C being the ordinary (one-sided) Laplace transform. The function A^q,(x) is given by 

OO OO OO 

MM = Yl ^\^n-a\{x) = e-^^/^' ^(-1)" J2.+^(2a;) + e-^^/^' ^^(-1)" J2„+„(2x) (22) 

n=— OO n=0 n=0 

where 7 = 2 — a. Substituting fl22l) in fl20|) . one gets that the Sumudu transform of Ma 
becomes 

[A^„(X)]5 = e-^^^ 5^(-l)" / e-'j2n+^{lXt)dt + 6""^^ ^(-1)" / J2n+c.{lXt)dt . 

n=0 n=0 

The use of the integral properties of the Bessel functions of the first kind [l^ (see pg. 386) 
gives the following expression: 

2n+a 



[Ma{x)], 



e '2' 



X' 



n=0 



Hi 00 



IX 



1 ^TTT^^r 

' * n=0 



VI -a;2 - 1 



IX 



1 i/(i-yr^)V Ai-v^ 



x^ 



X 



x^ 



X 



1/1 1 

+ 



2 vi-^' yr^ 



x^ 



'1 - Vl^ 



1 - 



(23) 
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By sending x — )■ 1/x in (123|) . and applying (12T!) . we obtain 



Ma - 



X 



\ ( X X 

+ 



2 V^' - 1 ^f: 



x^ 



=^ — Vx^ — 1 j + — Vx^ — 1 j 



and 



(x - V^^^)' " + (x - V^^^) 

2rx2 - 1) 



>C[A^„](x) 



(24) 



Hence the inter-particle statistical potential admits the following form, for many purposes 
easier to handle than (fTOl) . since it does not contain infinite sums: 

(25) 



(x — Va;^ — l)^ " + (x — \/ x'^ — l)' 



x^ - 1 



The correct potentials for the bosonic and fermionic cases are straightforwardly reproduced: 
using known results for the Laplace transforms 26|], we get 

2x 




e ^ 2 cosh X = 1 + e ; 



e ^ 2 sinh x = 1 — e 



(26) 



(27) 



Eq.( l25|) gives a closed formula for the potential in the case of semions [a = 1/2, or a = 3/2): 

' {x - v^^Y^^ + {x - v^^y'^^' 



where erf is the error function 



X2 - 1 



erf(V2x) (28) 



27|. 



Eq. (p5|) allows us to recover the second virial coefficient of a gas made of identical Abelian 



a— anyons, which is given by 



16|: 



1 



B2{a, T) = -A| (-1 + 4a - 2a^) . 



(29) 



The link between the 2"''^-virial coefficient and the statistical potential can be expressed in 
the form 

B,{a,T) = ^ dx [1 - e-^^(^)] . (30) 

Using Eqs.f l26|) -f l27|) -fl28 |) . for the three special cases a = 0, 1 and 1/2 (corresponding 
respectively to bosons, fermions and semions) one immediately finds 



x 2 roo 1 



X2 ■ 



(31) 
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B,{a = hT) = fj^ dx [1 - (1 - = +- 



a: 



(32) 



i?2(« = -,T) 



■\2 /-oo \2 poo -1 

^ rfx [1 - erf ( v^)] = ^ dyy[l- erf = - (33) 



as it should be. 

Finally, the effective 2-body statistical potential written as in Eq. f l2S]) allows us to easily 
recover the expression of the 2"'^-virial coefficient even for a general statistical parameter a. 
Indeed, by virtue of (125|) and the dominated convergence theorem, one has 

S2(«,r) ^_ 

" '0 



1 

2 y^ [1 - 



1 /"^ 

- lim / dx 

2 e^O 



- lim 

2 £^0 



(x — Va^^ — l) 
x2-l 



-1 



(a; — Va^^ — l) 
x2 - 1 



|a;=l+£ 



£-1 



(x — Va^^ — l) 



x^ 



£-1 



(x — Va^^ — l) 



- lim 

2 e-s>0 



v2 



2U-a 



+ (! + £- V2 



e + e 



2\a-l 



2e + e2 



^(-l+4a-2a2) , (34) 



that is just (1291). 

As a byproduct of Eqs. flT2l) . flTSl) . fl30|) . flM|l . we find an interesting integral property 
relevant to the bivariate Lommel functions (new, at the best of our knowledge): 



/ dx [l — 2 e^^ [i'" Ua{ix, ix) — U^-Jyix, -ix) ] } = — - + 2a 
Jo 2 



a 



(35) 



III. STATISTICAL POTENTIAL FOR NON-ABELIAN ANYONS 



In this Section we discuss the statistical interparticle potential for a two-dimensional 
system of S\J (2) NACS spinless particles. The NACS particles are pointlike sources mutually 



interacting via a non-Abelian gauge field attached to them [28|]. As a consequence of their 
interaction, equivalent to a non-Abelian statistical interaction for a system of bosons, they 
are endowed with fractional spins and obey braid statistics as non-Abelian anyons. 

Let's briefly introduce the NACS quantum mechanics 3Ch34|. The Hamiltonian de- 
scribing the dynamics of the A^-body system of free NACS particles can be derived by a 
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Lagrangian with a Chern-Simons term and a matter field coupled with the Chern-Simons 
gauge term js^: the resulting Hamiltonian reads 

N 

M,, 



i=l 



where Mj is the mass of the i-th particles, Vg. = ^ and 

" ail + 2^ 5 ^^"^^"i;^ ■ 

In formula fl36!) z = 1,...,A^ labels the particles, {xi,yi) = {zi + Zi, —i{zi — Zi))/2 are 
their spatial coordinates, and Q^'s are the isovector operators which can be represented 



by some generators in a representation of isospin I |32|. The quantum number I labels 
the irreducible representations of the group of the rotations induced by the coupling of the 
NAGS particle matter field with the non-Abelian gauge field: as a consequence, the values 
of / are of course quantized and vary over all the integer and the half-integer numbers, with 
1 = 1/2 being the smallest possible non-trivial value (/ = corresponds to a system of 
free bosons). As usual, a basis of isospin eigenstates can be labeled by / and the magnetic 
quantum number m (varying in the range — /, — / + 1, ■ ■ ■ , / — 1, /). 

Hence the statistical potential depends in general on the value of the isospin quantum 
number / and on the coupling k (and of course on the distance r and the temperature 
T). The quantity k, present in the covariant derivative is a parameter of the theory. The 



301 ]. In the following 



condition Auk = integer has to be satisfied for consistency reasons 
we denote for simplicity by k the integer Attk. 

For non-Abelian anyons, in analogy with iQ, the effective statistical potential can be 
related to the relative PF according to the following expression: 

Z^(«:,/,T)-Z^("-*-)(/,T) = ^ I d'pe-^^'/'^ j ciV [exp[-/3i;(K, /, r)] - exp[-/3t;(-^-)(/, r)]] , 

(37) 

where f (""*•)(/, r) refers to the system of particles with isospin / and without statistical 
interaction (k — )■ oo). The potential ?;("■*■)(/, r) can be expressed in terms of the potentials 
Va=Q{T) and Va=i{r) for the free Bose and Fermi systems (endowed with the chosen wave- 
function boundary conditions). 2'2(k, /, T) — '•""*'''(/, T) is the (convergent) variation of the 
divergent PF for the two-body relative Hamiltonian, between the interacting case in exam 
and the non-interacting limit k — )■ oo. 
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For hard-core boundary conditions on the relative two-anyonic vectorial wave-function, 
the quantity v^""^'\l, r) which enters Eq.f l37|) is given by the projection onto the bosonic/fermionic 
basis: 



1 



{21 



j=0 



1 + (-lp+2' 



-I3va=i{r) 



(38) 

in analogy with the procedure shown in 18N20| for the computation of the 2"'^-virial coeffi- 
cient. By using the results e~^'""=°^^^ = 1 -|- e~^^, e~^^°'^^^^'^ = 1 — e~^^, one then obtains 



exp[-^t;("-*-)(Z,r)] = 1 + 



-2x 



2/ + 1 



(39) 



Notice that this non-interacting quantity exactly corresponds to (—1)^^ times the statistical 
potential for a system of identical (2/)-spin ordinary particles (fulfilling the spin-statistics 
constraint) at the same temperature, similarly to what argued in 20| about the 2"'^-virial 
coefficient for the same system. 

In the interacting case (i.e. finite k), we can express the statistical potential in terms of 
statistical potentials of Abelian anyons: 



21 

e-^''('^'^'^) = .:r4TT^$:(2j + l^ 
^ 3=0 



{21 



e y 



(40) 



where Uj = [j{j + 1) — 2l{l + l)]/k, and v^.{r), v^^{r) are the potentials for the Abelian 
cuj-anyon gases respectively in the bosonic and fermionic bases, given by 



2e-^M^Ax) 



(41) 



and 



2e-^M^^+i{x) 



(42) 



Both (14T!) and (142!) are periodic quantities under the shift uj Uj + 2; it follows 



-l3v{k,l,r,T) 



2e- 



21 



(2/ + 1)- 



E + 1) 



j=0 



.1 V+2« 1 _ (_]\j+2l 

— M^^{x) + ^—^ — A^.,+i(x) 



• (43) 



Eq.( l43|l gives the statistical potential for a gas of NACS particles. 
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A. 2 -virial coefficient 



An useful application (and check, at the same time) of Eq. fH3|) consists in computing the 
second virial coefficient. The analogous of ( 130|) reads 



B2{kJ,T) 



\ 2 roo 



dx [1 — e" 



-I3v{k,l,r) 



(44) 



Substituting in its integrand both Eg. ( 1431) and the following decomposition of the unity 



(2/ + 1) 
one obtains for B2{k, I, T) 

21 



1 



3=0 



(45) 



2(2/ + 



roc 

— 5^(2j + l)/ dx 



\^ (1 - 2e-^M^Ax)) + "—^ (1 - 2e-^M^^+,{x)) 



By virtue of fl35l) one has then 

2/ 



j=0 



-7 + 7. -o7 + 



(46) 



-7 + ^i-7^.j 



where 7j = Ujmod2 and = {uj + l)mod2. This result matches with previous results 
reported in literature [19|, l20|, see in particular Eqs.(38) and (41) of [20 1. 



B. Minimum points 

Using Eg. (143 p . we can study if the gas of NACS has a "bosonic-like" , "quasi-bosonic" 
or "quasi-fermionic" behaviour, according to its characterization in terms of the statistical 
potential. In correspondence with the analysis carried out in Subsection III At we can ad- 
dress the problem of determining which points of the discrete parameter space {k, 1} are 
associated to the presence of an (interior) minimum point rcrit{k, I, T) for the statistical po- 
tential v{k,l,r,T) (which will be referred to as "bosonic" ones), which points correspond to 
a monotonically increasing v{k,l,r,T) ("bosonic-like" points) and which ones instead cor- 
respond to a v{k,l,r,T) monotonically decreasing in r (the "quasi-fermionc" ones). To this 
aim, let's exploit the limit behaviors of exp[—l3v{k,l,r,T)] for small distance {x <^ 1) and 
large distance {x ^ 1), which straightforwardly arise from Eqs. ([T5l) and (|T8|) . Notice that 
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at X = one has 

6-^^^''^'^='^ = j^^^ g (2j + l) [(1 + (-1)^-+^') 5(7„ 0) + (1 - (-1)^-+^') 5(7„ 1)]] , (47) 
where S denotes the Kronecker delta function and 7j = Ujmod2. For large distance it is 



t{j,k,l), otherwise. 



where 



and 

21 



21 

s{j, k,l)^Yl + 1)(-1)^+'' cos {u.n) 

j=0 



t{j, k,l)^Yl (2J + 1) sin(7,vr) [(1 + (-1)^-+^') (7, - 1) - (1 - (-1)^+^') {v, - 1)]] 

j=0 

with 7j = Uj mod 2, rjj = (uj + 1) mod 2. 

Our results can be summarized in the "phase diagram" shown in FigHJ in which we 
distinguish pairs of parameters {k, I) for which v{k, I, r, T) has an interior minimum point (in 
black), pairs for which the statistical potential is monotonically increasing in r (in magenta), 
and the remaining pairs (which are left blank). In this way the black points denote a "quasi- 
bosonic" behaviour and the magenta ones denote a "bosonic-like" behaviour, according to 
the classification operated in Section II to extract information from the statistical potential 
for Abelian anyons: one sees that for non-Abelian anyons one has mixed regions in which 
quasi-bosonic and quasi-fermionic behavior alternate, separated by regions dominated by a 
quasi-bosonic behavior. Bosonic-like behaviour occurs only for pairs {k, I) having one of the 
forms: {k generic, / = 0), {k = 1,1 integer), or {k = 2,1 even). 



IV. CONCLUSIONS 



In this paper we have studied the two-body effective statistical potential (which models, 
in the dilute regime, the dominant term of the statistical interaction between the particles) of 
ideal systems of Abelian and non-Abelian anyons, described within the picture of fiux-charge 
composites. In both cases we have derived closed expression for the statistical potential (in 
terms of known special functions) and we have studied its behavior. Asymptotic expansions 
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FIG. 4. Phase diagram in the parameter space {k, 1} for the family of non-Abehan anyon models 
under consideration. Points {k, I) in the black regions correspond to the presence of an interior 
minimum point for the statistical potential v{k, I, r, T) and are associated to quasi-bosonic behavior. 
Points in magenta correspond to a bosonic-like, monotonically increasing, behavior of the statistical 
potential. Finally, the remaining points are left blank and correspond to quasi-fermionic behavior. 

have been provided, and the second virial coefficients of both systems have been found using 
a compact expression of the statistical potential in terms of Laplace transforms. A phase 
diagram for the non-Abelian gas as a function of the Chern-Simons coupling and the isospin 
quantum number has been derived. 

In our study we have considered hard-core boundary conditions for the relative anyonic 
wave-functions: however, it would be interesting to use the results obtained so far in order 
to analyze the statistical potential in the more general soft-core conditions (both for the 
Abelian and non-Abelian ideal anyon gases). 

Acknowledgements: We wish to thank Andrea De Luca and Gianfausto Dell'Antonio for 
very valuable suggestions and discussions. 
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Appendix A: Definition of the used Bessel functions 

In the main text we used the Bessel functions of the first kind Jq, and the modified Bessel 
function of the first kind Jq,: their definition is respectively given by 

^ ^m\r(m + a + l) \2) 

m=0 ^ ^ 



and 



, , -(-) . (Al) 

ml r m + a + 1 V2/ ^ ' 

m=0 ^ ' 

The Lommel functions of two variables are defined in Eq. (5) of pg. 537 of [lj| and read 



m=0 



Vn{w, ^) = COS ( ^ + ^ + ^ ) + U2-u{w, Z) . 



Appendix B: Properties of the statistical potential 

In this Appendix we provide details and further informations on the statistical potential 
Va for Abelian anyons. 

We first give a derivation of Eq.f llOp : the relative PF ([7]) is given by 



2 l=-oo 



(Bl) 



and it can be rewritten as 



l= — oa \ / 



as pointed out in [11|. That allows for its comparison with the PF (jH]) for classical systems 
associated to a generic potential v{r), whence the result ffTOj) . 
We also observe that 



-I3vc,{r) _ 2 I 



00 

n=0 n=0 



(B3) 



The expression that follows here below is a possible closed form for the statistical potential, 

n 

but at the cost of using an integral representation given in formula (7), pg. 652 of It 
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stands for any complex number /i in the vertical strip — 1 < Re /i < a: 



nx) 



n=0 



?i=0 



[/•ia: fix 
y - 1) rft - r y - 1) Ji_«_^(t) rft . (B4) 

A simpler integral representation, valid for a G (0, 2), can be produced by using the following 
property ([3], pg. 540) of the bivariate Lommel function U: 

U^y{w,z) = J^^i{zt) cos !^^w{l-t'^)'^ fdt, Re(i^) > 

together with expressions fIBSp . f|T3l) and flAll) . The resulting integral representation is 



-^fa(r') _ 2 



X e 



/ cosh 




Jo 


12^ '\ 



(B5) 



In the final part of this Appendix we provide the derivation of the integral representation 
( IT^ . To this end, we use the following representation 25[ for the modified Bessel function 
of the first kind: 



1 2 cose n^n sini^TT 
cos v6d9 



= - / e 



-2 cosht— l/t 



TT 



Jo TT JO 

Then the summation term in ffTOD is: 



cit, arg|;z| < -, Rei^ > . (B6) 



-^|2n-a|(^) — -^2n+2-a(2;) + /2n+a(2) 



n=0 



n=0 



1 



vr 



cos (2?7, + Q;)0 + COS (2 



n — a 



n=0 



n=l 



where 

oo oo 

/(t, a) = ^ e-(2"+")* sin (2n + a)7r + ^ g-^^"-")* sin (2n - a)n 

n=0 n=l 

for t 7^ and 

/(O, a) = lim /(t, a) = (1 — a) sin avr . 



(B7) 
sinh[(l — a)t] 



sm avr- 



sinht 
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The first addend of the last integral representation is 



1 r 



vr 



cos (2n + a)(f) + cos (2n — q) 



+ cos (2n + 



.n=0 



a) 



n=l 



1 /■'^ 



vr 



COS </< 



.n=0 



/ d(l)e'''°'^ V cos (2?i + a)0 = - / 



^^^ ,5(2<A) + <5(2</)-27r) 
2 



e'""^ 2^ 



— ( e + e cos avr 
2 ^ 



cos (2n + a)(?l) + 

+ 00 

Re (ei"<Ae2m0) 

=— oo 

(B8) 



As a result of ([T0|), p37|) and ([BS]), one has then 

/3«.W = l + e-2-cosa7r-2''''"''-- 



vr 



sinh [(1 - a)t\ ,„,ht 
sinht 



(B9) 



By direct inspection this result, notwithstanding the hypothesis of validity for ()B6p . is valid also 
for Q at the extremes of the interval [0, 2], so that the derivation of ()14p is completed. 
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